Hydration of hydrophobic solutes treated by the fundamental 

measure approach 

O 

O , G.N. Chuev and V.F. Sokolov 

Sh , 

Oh. We have developed a method to calculate the hydration of hydrophobic solutes by 



< 



(D 



the fundamental measure theory. This method allows us to carry out calculations 
of the density profile and the hydration energy for hydrophobic molecules. An addi- 



^S^ . tional benefit of the method is the possibility to calculate interaction forces between 

O 

solvated nanoparticles. On the basis of the designed method we calculate hydration 



of spherical solutes of various sizes from one angstrom up to several nanometers. We 
C^ ■ 
"^ ■ have applied method to evaluate the free energies, the enthalpies, and the entropies 

■^.^ , 

C^ , of hydrated rare gases and hydrocarbons. The obtained results are in agreement 

I ' with available experimental data and simulations. 
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Q ■ Hydrophobic interactions play an important role in stabilization of various biomacro- 

O ■ molecular complexes including nucleic acidSjproteins,and lipids,because these complexes 



SjProte: 
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c^ ' contain a large number of nonpolar groups jjj, |2]. Despite the long history of studies of 

hydrophobic interactions, the theoretical treatment of the nature of hydrophobic interac- 

C ' tions is still incomplete. The main problem of the interactions is a complicated multiscale 

O' 

O ' character of these effects which can reveal as microscopic changes of water structure near 

- ^ , small hydrophobic groups, as well as conformations and aggregation of biomacromolecules 

><■ n 

;h ' at mesoscopic scales up to several tens of angstroms [3|]. The Monte Carlo (MC) and molec- 

ular dynamics (MD) methods are most frequently used for modelling molecular interactions 
in solutions |j]. However, their application to the problem of hydrophobic interactions of 
macromolecules demands huge computing expenses and in some cases it is essentially limited 
because of the specified multiscale character of these interactions. New methods based on a 



statistical treatment have been actively developed in last decade 



. The most 



suitable method for the specified effects seems to be the density functional theory (DFT) 
J9|, |lO| . The basic purpose of the approach is to construct the free energy functional of the 
system which depends on density distribution of liquids particles and intermolecular inter- 
action potentials. Within the framework of this approach there are a lot of various models 
connected with a concrete choice of the density functional |ll|, ll2, ll^, ll^ ■ 

In our opinion, one of the most perspective DFT models for calculation of solvation phe- 
nomena is the fundamental measure theory (FMT) |l2l . ll3L llSl Il6| which determines the free 



energy functional as the sum of the weighted contributions dependent on geometrical char- 
acteristics of fiuid particles. It automatically results in definition of the weighted functions 
which are responsible for the volume and the surface contributions to the solvation energy. 
This approach is intimately related to the scaled-particle theory (SPT) 121, llq] for the homo- 
geneous hard-sphere fiuid, and thus one expects that in the uniform limit it should reproduce 
the SPT results. The current status of this theory includes various generalizations of the 
scheme to different inter-particle interactions, binary mixtures and polydisperse systems, ap- 
plication to interfaces, wetting, confined geometries, porous media, and dynamical problems 
as well (see, review [15]). In this work we will calculate the radial distribution functions for 
rare gases and hydrocarbons, and thermodynamic parameters of their hydrophobic solvation 
on the basis of the modified FMT which treats the pressure of the liquid system correctly 
as distinct from the original FMT. In particular, we will determine the size dependence of 
the solvation energy and calculate interaction forces between two solutes within the limit of 
low concentration of the dissolved particles. The rest of this paper is organized as follows. 
The theory is described in section II, the obtained results and the concluding remarks are 
presented in section III. 

II. FUNDAMENTAL MEASURE THEORY 

The DFT is based on the unequivocal dependence between equilibrium density distribu- 



tion neg(r) and external potential Mextli") acting on a system 



equii 



The free energy of the 



system F[n] is related with the thermodynamic {Q) and the chemical (/i) potentials: 

f2[n, m] = F[n] — / drn(r)u(r), (1) 



where u{r) = fi — Uexti^^)- The equihbrium density neg(r) is determined by the minimuni 
Q[u] = Q[pi^,u] for the given temperature T: 

5n[n,v] _ 5F[n] 



= w(r). (2) 

neq{r) 



Sn{r 

In turn chemical potential fi is determined from boundary conditions, i.e., the equilibrium 
density should tend to the average density of homogeneous liquid neg{r -^ oo) -^ pj,. Thus, 
having the information about functional f2[n, m] and the method for evaluation of n(r), we 
can calculate all necessary equilibrium characteristics of the system. The density functional 
theory [9| solves the problem by searching the free energy F[n\, which consists of two con- 
tributions: the ideal {Fid[n]) and the excess (Fe^^fn]) free energies: 

F[n] = Fid[n] + Fe^n], PFid[n] = J n{r) \n[n{r)A' - l]dr, (3) 

where A is the de Broglie wave length, and (3 = [ksT)^^ is the inverse temperature. This 
free energy is related with the grand canonical functional, minimization of the functional 
leads to the equilibrium density: 

ne.(r) = p,exp[-/3.,.,(r) + ^-^ - ^^^^^ ]• (4) 

Thus, if the functional -Fe2:[n] is known we can calculate the density profile ?T,(r) and then all 
required characteristics of the solvation. 

There are many ways for constructing Fei:[n], most of them use the data on functional 
derivatives dFe;j./dn{r) and d'^Ff.x/dn{r)dn{r'). For Lennard- Jones (LJ) fluids the excess 
free energy is decomposed into the contribution from a reference system of hard spheres, 
and the free energy due to attractive interactions 

Fe.[n]=F,,M+F,,Jn]. (5) 

The attractive interactions are usually treated by the flrst-order or the second-order pertur- 
bation theories |20,|2l|. In the mean fleld approximation the attraction potential Uatt{f) is 
considered as a perturbation which gives the contribution to the free energy 

Fatt[n] = 2J j^^^^^~ Pb\Uatt{r - r')[n(r') - p^]dv'dv. (6) 



Various functional for inhomogeneous hard-sphere (HS) fluids have been developed. One of 
such methods is the FMT [l^ |l3|, |l5|, |l6| in which the excess free energy is calculated by the 



use of coarse-grained or smoothed densities. We note that originally this method has been 

ormulated for HS liquids. But later this method has been applied to spheroids of rotation 

"1 

22| and also for various repulsive and attractive potentials 

repulsive contribution to the excess free energy is written as 



2J,|2^. In HS liquids the 



f3Fhs[n] 



$[ni(r)]dr, 



(7) 



where variables nj(r) are determined as weights 

ni{Y) = f dr'n{r') w^'\r-r'), (8) 

of the density n(r') averaged with weight factors w^^\v — r') depending on the fundamental 
geometrical measures of fluid particles, such as volume, surface, etc. The original Rosenfeld 
formulation 12, ll^ utilizes the following weighting functions 

w;(3)(r) = e((T/2-r), w^^\r)=6{a/2-r), w^^\r) = w^^\r)/{2TTa), (9) 



r 
where S{r) and G(r) are the Dirac delta-function and the Heaviside function, respectively, 
a is the diameter of a solvent particle. These weight factors determine weighted densities 



no(r), n„i(r), n^^r), and latter is t 



Within the FMT framework 
weight densities raj(r): 



12 



le local factor of packing. 
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the function $[nj] is determined through six 



$k; 



-"° ^"(' - "^) + 1-^3 + 2Ml-n3)^ • ('°) 

Using these equations and (jlj) we obtain the equilibrium density 

negir) = Phexp[-(3uextir) + Y.{J T—F^M^ - r')dr' ^^— ^)]. (11) 

Although the original FMT and the SPT is restricted by the HS fluids, there are no limi- 
tations for solute-solvent potential Mea:t(^)- Hence we may apply the FMT to treat the LJ 
solutes. The information on the weighted densities also allows us to calculate the mean force 
potential W{r) determining the interaction force between two solutes in an infinitely diluted 
solution 



W{r) 



Kn(r) 



E< 






m 



where Kiu(r) is the direct intermolecular interaction potential. In turn the excess part of 
thermodynamic potential Qex determining the solvation energy is calculated as 

i ^ 

Since experiments are most commonly done at fixed pressure p, it is convenient to introduce 
the decomposition of the hydration chemical potential into the excess solvation entropy 
AS* and the excess solvation enthalpy AH achieved by the use of an isobaric temperature 
derivative J26|, 

AS = f^^) , AH = A/.,, + TAS. (14) 

III. RESULTS AND DISCUSSION 

A. Application to bulk water 

The FMT reduces to the SPT in the limiting case of homogeneous liquid, determining 
the pressure p and the surface tension j^ of the liquid at a planar wall ^] as 

(5<l>(r -> oo) ^ (5$(r ^ oo) ,^ ^, 

Oris dn2 

In the general case the pressure phs is too high to describe fiuids under normal conditions, 
for example, it yields 8000 atm for the effective HS diameter a = cr^„ = 2.77A corresponding 
to water at 25°C |28l |29|. On the other hand, formula (fT3|) underestimates the surface 
tension 7^^ with respect to the experimental values. To exclude these drawbacks both the 
thermodynamic parameters are considered as the fitting ones in the modified SPT models 



(see, for example J30|,|3l|). Thus we are to modify the FMT to obtain the realistic estimate 
for the surface tension. There are a lot of ways to do it, we use the simplest one by exclusion 
small distances at r < Vcut in integrals ©, (0), and (jH)). As a result, the formulas for the 
thermodynamic parameters are also modified: 

5$(r = (X)) ,5$(r = rc„t) ^. ^^(r = 00) 6^{r = r^ut) .,„x 

PPhs = 7 7 , P7oo = 7 7 • (Id) 

oris ons dn2 on2 

We chose Vcut to fit ^y^ by the experimental value Je^p- Below we will indicate that the 
modified pressure phs also reduces significantly by three orders and plays a minor role in the 
hydration of hydrophobic solutes. 
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First, we employ the method to evaluate the thermodynamic parameters of bulk water. 
The HS diameter for water was chosen a^ = 2.77 A and density of water was the PtO"^ = 0.7. 
Such choice of the water diameter has been motivated by the solubility experiments of 
Pierotti ^28|. To estimate the attractive contribution (jH)) we have used the LJ parameters 
corresponding to the SPC/E model of water 221 • For the solution of the equation (fTTjl 
we used the Picard iterative algorithm. On every fc-step of iteration it was necessary to 
calculate weight densities nj(r), which are related with neq{r) through integrals (jHl). For 
improvement of convergence we used the algorithm based on mixing of parts of the previous 
and new iterations n^^^ = Anf„ + (1 — A)n^„4, where A is the mixing parameter dependent 
on the bulk density. The step of integration made O.Olo"^, and the number of points of 
integration is A^ = 2^^. These numerical parameters provide the relative precision of density 
profiles up to 10~^ and for the chemical potential up to 0.1 kcal/mole. Using relations ()13p 
with the cut-off radius Vcut "we have evaluated the excess chemical potential depending on 
the diameter of the HS solute: 

A/ig^(cr„) = fiez(0-„) - QexicTu = 0), (17) 

where cr^ is solute diameter. We used Vcut = 0.42(7,„ to fit the calculated surface tension to 



the experimental value •y^^p = 102 cal/(mole A ) [23|- Figure 1 represents the dependence of 
the excess chemical potential A/ig^(cr„)/47r(T^ obtained by this fitting. Apparently from the 
figure our calculations are very well compared with the MC results 3^ and the calculations 
received on the basis of the information theory .35] • The figure shows that on the site 
from zero up to 4 A the curve behaves almost in the linear fashion. It indicates that at such 
solute sizes the surface effects do not yield the appreciable contribution to the excess chemical 
potential and the volume contribution has crucial importance. Whereas for the particles, 
whose radius is more than 10 A, the contribution of the volume component decreases and 
the surface component grows. 

We have also calculated the Tolman length 6, which is the surface thermodynamic property 
of the water vapor-liquid interface (the distance between the equimolar surface and the 
surface of tension). For this purpose we use the relation 

^ f R-25 \ 

Ap,, = 87ri?%^ L^_^^ j , (18) 

where R = au/2 is the solute radius. As a result, we have obtained 5 = 0.92 A which is 
agreed the MC simulations of the SPC/E water '^^. This value is a little bit more than 



that calculated in [30| {5 = 0.9 A). The pressure of bulk water is estimated as 

P = Phs + ^ j Uatt{r)dY. (19) 

Table 1 lists the data on the above thermodynamic parameters obtained by the FMT with 
and without cut of the integration range, as well as the data derived from MC simulations 
3J| . Although the pressure obtained by our procedure exceeds the simulated one by an order, 
it does not yield the significant effect on the dependence A/ig^(o"„), since the pressure effect is 
minor for water under normal conditions. These results hold, however, only at temperatures 
near 25° C, since the vapor-liquid interfacial tension of water decreases monotonically with 
temperature increasing, while the excess chemical potential of hard sphere solutes exhibits 
a maximum with temperature increasing. 

B. HS and LJ solutes 

On the basis of the designed method we have carried out calculations for various hy- 
drophobic objects. We have considered the two models of hydrophobic solvation. The first 
of them consists in that the solute is modeled as hard sphere. In framework of this model we 
have obtained the dependence of the excess chemical potential on the radius of HS solute. 
Figure 2a shows the comparison of the results received by FMT and by Monte Carlo simu- 
lation 29|. The discrepancy between two results are practically missing. The FMT results 
a bit underestimate the excess chemical potential for hard-sphere solute which radius more 
than 4 A. To justify our results we have also calculated the enthalpic and entropic contribu- 
tion into the chemical potential for the HS solute. Figure 2b shows the comparison of the 
results calculated by the modified FMT with MC results 29j. The main difference between 
the FMT and the MC results is that the FMT as the SPT yields a monotonic decrease of 
entropy versus the HS diameter, while the MC data indicate more complicated behavior. 
On the other hand, the magnitude of the excess chemical potential for the HS solutes which 
radius is a less 6 A received by MC a little bit exceed the FMT results. 

The LJ model consists in that interaction between solute and solvent molecules is realized 

by the LJ potential 

cr \12 /a 



u...{r)=Ae[[^-fL^ -(^r^j], (20) 

where Suv and auv are constructed from the corresponding parameters e„, (T„, e^ and a^ by 
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Lorentz-Berhlot mixing rules. Figure 3 shows the two examples of the radial distribution 
functions (methane and neon) calculated by the modified FMT and derived from the sim- 
ulations J2IJ, 1^ ■ The discrepancies between the first maxima of the peaks do not exceed 
10%and the widths of the peaks are practically coincide. The reason of this difference un- 
derlies that the water molecules are modeled as hard spheres and the peaks become higher 
and in many instances a bit narrower. Since we have modeled the solutes as spheres this 
model seems to be more essential for calculation the excess chemical potential of the atomic 
solutes. Besides we will apply this approximation to treat the hydration of hydrocarbons 
also. 

C. Inert gases and hydrocarbons 

For the evidence of the efficiency of the FMT we have calculated radial distribution func- 
tions for linear, branched, cyclic hydrocarbons (methane, ethane, butane and et al.) and for 
rare gases in water. The LJ parameters are presented in Table 2. On the basis of HS model 
we have calculated radial distribution function and evaluated the excess chemical potential 
for the hydrocarbons and rare gases. The calculation of the excess chemical potential has 
been carried out with the use of the perturbation theory (PT) to take into account the at- 
tractive contribution of the solute-solvent interactions. In this case we have calculated it by 
using equation ©• The attractive solute-solvent contribution to the hydration free energy 
has been estimated as 

Fuv = Pb ghs{r)Uatt{r)dr, (21) 

w.e. .,.) . co.e,a«on .n..o„ of „a.e. .o.ec.e. _M .he HS .o,u. Q. Fo. .M. 
contribution we have used the Weeks- Chandler- Anderson decomposition for the LJ potential 



into the attractive and the repulsive parts [2 

Uatt{r < 2^/V„„) = -e, Uatt{r > 2^/V„„) = 4e, 









(22) 



where auv and Euv are the LJ diameter and well depth, respectively. We denote this approx- 
imation as the LJPT model. 

Note that the attractive contribution of solvent-solvent interactions energy decreases from 
15 to 1 percent with increasing the solute radius from 0.5 to 20 A. Opposite, the attractive 
part of solute-solvent energy increase with increasing the solute radius. Its magnitude is 
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about 80% of the excess chemical potentiaL Table 2 shows the calculated values of the 
excess chemical potential, the enthalpies and the entropies for hydrocarbons. The magnitude 
of the excess chemical potentials for small hydrocarbons (methane, ethane and propane) 
is hardly different from experimental results. But the difference of the excess chemical 
potentials for large solutes becomes more significant. The main reason of it consist in that 
the solute-solvent attractive contribution for large solutes are essential and it's necessary 
to know the realistic value of Eu- Since hydrocarbons are molecular solutes, unfortunately 
we can't apply an unique LJ parameter e^v correctly. In much the same way we have 
calculated the excess chemical potential for rare gases (He, Ne, Ar, Kr and Xe). At first, we 
have calculated the energy of cavity formation and two corrections which take into account 
attractive contributions to the solute-solvent and solvent-solvent interactions. Figure 4 
shows the insignificant difference between the experimental values of the excess chemical 
potential and received by the modified FMT. The agreement between the experimental and 
the calculated data is strongly quantitative. Table 2 shows that the excess chemical potential 
with corrections of attractive contributions of the solvent-solute and solvent-solvent energy 
depends on solute radius in the "U" form J2S|- The left and right parts of sign "U" put 
together the rare gases and the hydrocarbons, respectively. 

We have also applied the LJ model and utilized the Lennard- Jones potential for solute- 
solvent interactions. The discrepancy between the experimental and the calculated values 
of the excess chemical potential for the rare gases and hydrocarbons is practically missing. 
Table 2 shows the LJ parameters of the rare gases and hydrocarbons and the results of the 
calculations. We have carried out the analysis of decomposition the excess chemical potential 
for hydrocarbons and rare gases on entropic and enthalpic parts. Using Eqn. ()14|) . we have 
derived the excess chemical potential to two parts Ai7 and —TAS. Table 2 lists the data 
on the calculated and the experimental excess chemical potential, the enthalpies, and the 
entropies of hydration of the hydrocarbons from methane to hexane, the branched hydro- 
carbons (2-methylpropane, 2-methylbutane, and neopentane), and the cyclic hydrocarbons 
(cyclopentane and cyclohexane). The FMT calculations by the LJPT and the LJ models 



are labeled as LJPT and LJ, respectively. We note similarly [3S| that the obvious feature of 
the entropic and enthalpic terms that its are greatly larger in absolute value than the excess 
chemical potential. The hydration enthalpies are large and auspicious, and the hydration 
entropies are large and unfavorable. The entropic terms are marginally larger in absolute 
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value than the enthalpic terms resulting in the unfavorable but small hydration free energies 
of the hydrocarbons. It is recognized that this behavior is typical of hydrophobic hydration. 
Solvation of apolar compounds in most other solvents, in fact, is usually accompanied by 
smaller enthalpic and entropic changes. Table 2 indicates the difference between the exper- 
imental j4i| and the calculated enthalpies and entropies. The discrepancy in enthalpy does 
not exceed the 7 kcal/mol for all hydrocarbons. The calculated excess chemical potential 
of methane is about of 3% more positive than the experimental value (2.01 kcal/mol). The 
discrepancy between the experimental and calculated hydration enthalpies for normal hy- 
drocarbons is less significant then the similar quantity for the higher hydrocarbons. The 
difference in the enthalpy and the entropy partially cancels each other resulting in a smaller 
discrepancy in the excess chemical potential. The calculated and experimental hydration 
free energies of the hydrocarbons are positive and practically do not increase with solute 
size. 

The fact that the excess chemical potential of ethane is less than methane is well repro- 
duced by our calculations. The modified FMT overestimates the value of enthalpy change 
and too highly underestimate entropy loss for molecules from methane up to ethane. Both ef- 
fects bring in less favorable hydration free energy of ethane. It appears, because the current 
model of hydrocarbons should take into account a larger benefit in favorable hydrocarbon- 
water interactions in going from methane up to ethane without the further loss of entropy. 
The hydration enthalpies of the normal and linear hydrocarbons have a tend to increase as 
solute size rises. The calculations reproduce qualitatively the effect for the LJPT model. 
The theory overestimates the magnitude of the hydration enthalpies. Fortunately, the dif- 
ference between calculated and the experimental results are rather closer for the LJ model. 
The analogue behavior is seen for hydration entropies. 

Similarly such decomposition the excess chemical potential into entropic and enthalpic 
part has been carried out for rare gases. Table 2 shows the discrepancy of solvation free 
energy for all rare gases to be negligible. The difference between the calculated and the 
experimental excess chemical potential for all atoms does not exceed 0.33 kcal/mol at tem- 
perature 298 K. The calculated solvation entropy are not so close to the experimental value, 
unlike the calculations of the solvation enthalpy which less overestimate the contribution of 
it to solvation free energy. The theory underestimates the values of solvation enthalpy and 
overestimates the solvation entropy for all rare gases. 
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The hydrophobic effect is frequently connected to characteristic temperature dependences 



41 



42J | . One of the most surprising observations is the entropy of transition convergence of 
nonpolar molecules from gas phase or nonpolar solvent into water at temperature about 400 
K. We have made the calculations to show that the FMT is able to predict the temperature 
convergence of entropy both qualitatively, and quantitatively correctly. The calculations 
have been carried out at several temperatures along the experimental saturation curve of 
water. Using ()14|) we have obtained the solvation entropy by taking the derivative of the 
chemical potential along the saturation curve. Figure 5 shows the temperature dependence of 
the entropy for the different solutes for two cases of calculation. In the first case (Fig. 5a) we 
calculated the excess chemical potential without the fact that the diameter of water decreases 
with temperature increasing. The entropies are large and negative at room temperature for 
all the solutes and decrease in magnitude with increasing temperature. The temperature 
dependence of entropies is approximately linear with slopes increasing with the increasing 
solute size. Moreover, the entropies converge at about 400 K to approximately zero entropy, 
although at closer inspection the temperature range of the convergence region is several 
10 K and the entropy is not exactly zero at convergence. In the second case, we have 
taken into account the dependence solvent diameter on temperature J43|. Figure 5b shows 
that the point of entropy convergence has shifted to region where temperature and entropy 
magnitude is about 470 K and -2.5 cal/(mol K), correspondingly. The convergence region 
has become a bit wider. It is significant that taking into account the contributions solute- 
solvent interactions has changed both the point of entropy convergence (about 500 K) and 
the width of the convergence region. 

D. The mean force potential for colloids 

We have to note the one more benefit of the FMT. The theory allows to calculate the 
depletion forces between two the solutes surrounded solvent particles. For macroscopic 
objects there are relations for calculation of depletion force for large solutes, depending 
on the distance between them. The solutes are located in environment of small solvent 
particles. Viewing hard spheres, Bradley and Hamaker have received the relations which 
take into account the pair interactions between macroparticles solvated in a simple fluid. 
Following this approach, the depletion potential W and the depletion force F between two 
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macroparticles have been calculated depending on the distance h between the solutes 



Q: 



.,„)._!lgl, «M.)^-^. (23, 



o 



2 , fs^-A' 
+ ^ + ln 






where s = 2R + h/R. We have calculated the depletion force F{h) on the FMT basis with 
the use of the correlation functions obtained before. Apparently from Fig. 6, the Hamaker 
and the Bradley approaches badly feature the behavior of depletion forces near to the solute, 
yielding strongly underestimated values, and not taking into account the oscillating character 
of it. The Hamaker and the Bradley approaches adequately predict the depletion force values 
only in asymptotic distances between the dissolved particles. Unlike them the FMT allows 
us to take into account attractive interactions between solutes a bit less than one nanometer. 
We have also compared the results of MD simulations 45| and that obtained by the FMT. 
Figure 6 shows that the first peak of FMT results is a bit narrower and higher than MD 
simulation but localization of curves' zeros almost coincides. The reason of narrowing of 
the peaks seems to be the modelling of water molecules as hard spheres not as particles 
interacting via Lennard- Jones potential. 

IV. SUMMARY 

In this work we have used the FMT for the quantitative description of the hydrophobic 
phenomena on the basis of the density functional theory. As a result, we have received 
profiles of radial distribution functions for isolated solutes in a hard-sphere fluid interacting 
with solute by Lennard- Jones potential. Using the distribution function profile, we have 
constructed the dependence of the excess chemical potential on the radius of the spherical 
solutes. To fit the properties of bulk water, namely its surface tension, we have modified the 
FMT by cutting the radius of integration at small distances. The excess chemical potential 
has been calculated for several systems, i.e., when interaction between the solute and solvent 
is simulated as the hard sphere or as the Lennard- Jones potential. We also have shown that 
the obtained distribution functions reproduce with good precision the oscillating behavior 
of depletion forces for the particles dissolved in fluid. We have applied method to evaluate 
the free energies, the enthalpies, and the entropies of hydrated rare gases and hydrocarbons. 
The obtained results are in agreement with available experimental data and simulations. 



13 

We conclude that the original FMT which rigorously should be applied to liquids where all 
the interactions are described by HS repulsive potentials is also well relevant for realistic 
water-solute potentials. However, such success requires the hard sphere water diameter to 
be used as an adjustable parameter. The recent extensions of the FMT to soft potentials 
indicate that the soft FMT is capable to predict the solvent structure for soft repulsive 



and attractive interaction potentials |46l . |47| . In this case the FMT yields a systematic 
way to generalize the treatment for hard bodies to soft interactions. Another bottleneck of 
the current implementation is the spherical shape of solutes, since the realistic applications 
should treat the three-dimensional solute structure. There are no restrictions to generalize 
the above approach to the three-dimensional case, the examples of such generalization are 
presented in (48L l49L l5fll |. however such calculations demand special algorithms and more 
sophisticated computations. Thus, we think that the FMT can provide a promising basis 
for the accurate study of hydrophobic molecular solutes. 
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Table 1. Thermodynamic parameters of bulk water obtained by the simulations and by 

the FMT. 





MC [34] 


FMT 


Modified FMT 


Pp<tI 


5.17x10"^ 


6.05 


4.94x10-3 


(^i< 


1.36 


1.06 


1.32 


5/ayj 


0.32 


1.59 


0.33 
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Table 2. LJ parameters of rare gases and hydrocarbons J38l . |39|, the experimental |40| 
and calculated data on of the excess chemical potential, the enthalpies, and the entropies of 
their hydration. 





A/ig^., kcal/mol 


-Aif , kcal/mol 


—TAS, kcal/mol 


Molecule 


o-«, A 


£u/kB 


exp 


LJPT 


LJ 


Exp 


LJPT 


LJ 


Exp 


LJPT 


LJ 


Helium 


2.63 


6.03 


2.75 


3.00 


2.76 


- 


0.49 


0.16 


- 


3.48 


2.91 


Neon 


2.79 


35.7 


2.67 


2.62 


2.62 


0.35 


1.23 


0.91 


3.02 


3.84 


3.53 


Argon 


3.41 


125 


2.00 


2.21 


1.91 


2.38 


3.09 


2.99 


4.38 


5.30 


4.89 


Krypton 


3.67 


169 


1.66 


1.98 


1.54 


3.20 


4.05 


4.12 


4.86 


6.04 


5.67 


Xenon 


3.96 


217 


1.33 


1.60 


1.38 


3.85 


5.22 


5.40 


5.18 


6.83 


6.79 


Methane 


3.70 


157 


2.01 


2.07 


2.11 


2.70 


3.97 


1.70 


4.71 


6.04 


3.80 


Ethane 


4.38 


236 


1.84 


1.61 


1.69 


3.90 
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Figure captions. 

Fig. 1. Dependence of the excess chemical potential/surface area of the solute on its 
HS radius. The solid line is plotted for the modified FMT, the dashed and dash-dotted 
lines for the modified SPC |3Cl| and the IT |l2^ , respectively, the triangles correspond to MC 
simulations (SPC/E water) |84| . 

Fig. 2. The excess chemical potential, the entropy, and the enthalpy calculated by the 
MC J29| and by the modified FMT: a) the excess chemical potential (solid and dashed lines 
plotted for the FMT and the MC results, respectively), b) entropies (black circles) and 
enthalpies (white circles). Symbols with lines and without of them denote the FMT and the 
MC results, respectively. 

Fig. 3. Solute-proximal water oxygen radial distribution functions for neon and methane. 
The solid line corresponds the modified FMT, while the dashed one to the MD results for 
Ne y and the MC data for Me |3o|, respectively. 

Fig. 4. The calculated and the experimental j^ (solid and dashed lines, respectively) 
data on the excess chemical potentials (triangles), the enthalpies (black circles) and the 
entropies (white circles) for the rare gases. 

Fig. 5. Hydration entropy of HS solutes with radii corresponding hydrocarbons (see 
legend) as a function of temperature along the saturation curve of water in the case when 
the water diameter is independent on temperature (a), and temperature dependent (b) like 
as in y. 

Fig. 6. Dependence of the depletion force on the distance between two solutes with 
radii Ru = 8.35 A and (3euv = 5, /3ev=l, a^ = 3.2 A. The sohd line corresponds to the 
FMT results, the dashed one to the MD simulations |,.45,j, while triangles and circles to the 
calculations by the Hamaker and the Bradley formulas. 
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